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Abstract 



The Modified Newtonian Dynamics (MOND) is presented here, as well as a theory that 
can be linked to it: the theory of the Aether, a four-vector field breaking Lorentz invariance. 
The form of its Lagrangian is studied, then basic equations of the theory are rederived in a 
detailed way, and calculated for different metrics, exploring the impact of non-zero spatial 
terms of the Aether. A brief attempt of making the Aether Lagrangian depend on a scalar field 
is presented. An analytic solving of a galaxy model with an external field is described, which 
highhghts the MONDian external field effect that breaks the strong equivalence principle. 



Resume 



On presente ici la gravitation Newtonienne modifiee (MOND) puis une theorie qui pent 
lui etre reliee, celle de I'Ether, quadrivecteur brisant I'invariance de Lorentz. La forme du 
Lagrangicn consideree pour I'Ethcr est etudiee, puis Ton rederive en detail Ics equations de 
base de la theorie, et on les calcule pour differentes metriqucs, cn cxplorant I'impact de 
termes spatiaux non nuls pour I'Ether. On presente ensuite unc tentative de faire dependre 
d'un champ scalaire le Lagrangien de I'Ether. On s'interesse egalement a une description 
analytique d'un modele de galaxie sous un champ exterieur en gravitation Newtonienne 
modifiee, faisant ainsi ressortir I'effet de champ exterieur de la MOND, qui brise le principe 
d'equivalence fort. 
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Introduction 



Gravitationnal physics is a great field of research and questionning for physicists. Newton's 
second law and General Relativity are no more sufficient today to account for astronomical 
observations without having to resort to uncertain forms of matter. Astronomical observa- 
tions show that the mass observed in the universe is far too small to account, for instance, 
for accelerations of objects in the outskirts of galaxies or for flat rotation curves of galaxies 
beginning at a certain radius. A possibility to solve this problem is to resort to Dark Matter, 
named DM hereafter, which would consist of particles whose nature we can but speculate 
about today, except that they of course do not couple to photons. By adding DM, we can 
"artiflcially" raise gravitationnal fields, and by a clever fitting, adding to galaxies a halo of 
DM whose density decreases in we can get flat rotation curves which are observed. 

We can also think in a different way: maybe we are looking for something that does not 
exist. DM is indeed needed for observations to fit with our current theory of gravitation, 
but maybe this theory is not always correct. Newton's second law and General Relativity, 
that is build to reduce to it in a limit of low speed and weak slowly variating potentials, are 
successful in many cases. They describe very well our Solar System and phenomena in it 
such as the precession of the perihelion of Mercury or the Shapiro effect. But maybe it is not 
always true for all scales. 

I worked during this internship under the supervision of Dr HongSheng Zhao working in 
the Physics and Astronomy department at the University of St- Andrews in Scotland. He's 
interested in DM and various aspects of the MOdified Newtonian Dynamics: galactic dy- 
namics related issues, gravitationnal lensing, etc... He and his two Phd students Xufen Wu 
and Garry W. Angus form the MOND group at the University of St-Andrews. I had no 
predefined goal at the beginning: my supervisor wanted me first to get familiar with Aether 
theories poping up in the relativistic side of MOND and to rederive equations of the theory so 
as to carry out further work. As a member of this group, I read quite a lot about MOND and 
got in touch with some of the related work, and even worked on a specific MOND problem 
(Chapter 2). Meanwhile, I mostly focused on the Aether which is by itself linked to a huge 
part of current theoretical physics. 

The next section is a short presentation of MOND, the first chapter is a study of Aether 
theories that can be linked to MOND and a report of my work about it, and the second 
chapter presents a "real MOND problem" linked to galactic dynamics and the "External 
Field Effect", which is a very important and controversial aspect of MOND, because it 
breaks the strong equivalence principle. 

1 MOND 

The Modified Newtonian Dynamics emerged in the eighties. Milgrom wrote in 1983 several 
papers (see [H], [E] and [K]) to introduce it as an alternative to dark matter. The problems 
encountered seem to occur for low accelerations, so Milgrom proposed to modify the Newton's 
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second law this way. 

rUgaii ( — ) = F (1) 

The value of ao, about 1 A was estimated by Milgrom in different ways, amongst 
which the fitting of rotation curves. 

The function fi appearing here is the MOND interpolating function, which conveys the 
transition between the Newtonian and MONDian regimes. The exact form of n is not clear, 
different attempts have been made to fit observationnal data, but its asymptotic behaviour 
is: 

• fi{x) = 1 when a; 1 such that Newtonian gravity is recovered as expected when the 
gravitationnal field is far bigger than the acceleration scale oq. 

• fi{x) = X when x <^ 1 such that the formula will reduce to gg = ctogN- 

But it is not clear whether just gravity or inertia must be changed. We can choose to 
change just the gravitationnal field this way: 

gf. [I) = gN (2) 

Is the whole physics in low acceleration modified? Some recent papers (see [S] and [7]) 
suggest that tests could be performed in the following years to try to check this on the Earth. 
This would anyway not be easy because it requires a high precision. It could rule out the 
MOND modification of inertia. 



Such formulas are empirical. They describe lots of observations very well, but do not 
consist a theory. Bekenstein and Milgrom wrote in 1984 p] an action for non-relativistic 
MOND, so that by variating it, one can find a modified Newton-Poisson equation. One 
can derive the Poisson equation by differenciation with respect to (w.r.t. hereafter) the 
gravitationnal potential of the following action: 

Sn = - j d^r [p^N + {SnCy^ (V$7v)^] V^^tv = 47rGp (3) 
whereas the following modified action: 

S = - [ rfV[p$ + (SttG)-' alF [ 

will generate the modified Poisson equation: 



(4) 



IV<I>I 

V.(/i(J ^)V$) = A-kGp, fi{x) = F\x^) (5) 

ao 

With this action, momentum, angular momentum and energy are conserved (see [2j). The 
function p and the constant ao are however still put ad hoc, with no physical basis. 



4 



2 MODIFICATIONS OF THE EINSTEIN EQUATIONS 
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2 Modifications of the Einstein equations 
2.1 Relativistic side of MOND 

A lot of work has been done to build a covariant theory around MOND. This theory would 
better be derivable from an action, namely to get conservation laws. Generally speaking, one 
builds a function of the dynamical variables of a theory, an action, that when variated w.r.t. 
these dynamicals variables gives the dynamic equations. Such a formalism is interesting be- 
cause conservation laws follow then directly from symetry principles. 



The Einstein-Hilbert action used for GR is: 

R 



S = d XV— ^ 



+ Sm (6) 



where: 



• g is the determinant of the metric gaf3- This metric defines the geometry of spacetime. 
For two events of spacetime Pi = (x") and P2 = (xg), the infinitesimal interval ds 
between them is given by ds"^ = g^pdz'^dx'^^*"' , where the dx" are the coordinates differ- 
ence between the events: xf = (x°) + dx"". The signature I took for all calculations in 
this report is (— , -|-, -f, +). 

• i? is the Ricci scalar. There are different conventions for the Rieman tensor R'^ ^au- I 
took for all the report: 

( V.V, - V,V.) = R^ (7) 
R fiau = d^T^j^ — duT^^ + r^j^r^^ — V^^Tj^p (8) 

The Ricci tensor R^,^ and the Ricci scalar R are:: 

RtMu = R^ ^iXu, R = g'^'^RtJ.iy (9) 

These tensors describes the curvature of spacetime. Q 

• Sm is the matter action that describes the matter distribution. 
The light speed is c = 1 in all the report. 

Differenciation of this action with respect to the metric gives the Einstein equations: 

1^ = ^ G., = SvrGTX*^'- (10) 

with : 

• Gai3 = Rai3 — \ga(3R'- the Eiusteiu tensor 



^The covariant derivatives and the Christoffel symbols are defined the usual way: for a tensor T^'" , 

^,T^ - = d,T^ - + r(;^r( - + ... - r^^r^' - - with rj;^ = {d^g^p + dpgx. - d^g^p) 



2 MODIFICATIONS OF THE EINSTEIN EQUATIONS 



6 



• T^^**^'^: the stress-energy tensor of matter defined by: 5Sm = ~l I d'^Xy/^T^y (x) Sg'^^ (x) 
This tensor describes the matter distribution. 

These Einstein equations hnk the geometry of spacetime (by way of the curvature of its 
metric) to the matter distribution (by way of its stress energy tensor). 

If we consider terms depending only on the curvature of space-time by the intermediate 
of the Ricci tensor as in GR, it can be shown that the Newtonian gravity will always be 
recovered [1], so we can't hope to recover MOND that way. We must add new degrees of 
freedom. After former attempts to build a fully covariant theory which encountered problems 
such as violation of causality, in 2004, Bekenstein proposed a theory named TeVeS, for Tensor 
Vector Scalar (see P) in which he considers three distinct degrees of freedom: a tensor (the 
metric), a vector and a scalar field. Zlosnik et al. noticed in 2005 that this theory could be 
reduced to a purely vector-tensor theory (see |23]). 

I did not work on TeVeS at a whole, but just on the Aether, the four-vector field that 
appears in it. The use of such a vector is embedded in a series of attempts to build Vector- 
Metric theories of gravity. Such a vector was considered by Will and Nordvedt in 1972 [22] . 
amongst other means of exploring the possibilty and the impact of "preferred frames" in 
gravitationnal theories. Once again, this upsets physics because it breaks Lorentz invariance, 
which reads that there is no preferred frame. Lorentz invariance can't be tested uniformly 
because a parameter of the Lorentz group is unbounded, so we don't know if it holds at 
each scale. What's more, this invariance leads to divergences in quantum fields theories, so 
exploring breakings of it is "allowed" and tempting. 

Breaking the invariance globally by choosing for instance a universal preferred rest frame 
or fixed background tensors is not satisfying because one can't this way preserve general 
covariance that is required for the Einstein equations to hold. One should give up GR or 
any modification of it, which is not theoretically appeahng (see [9j). One therefore tries to 
break this invariance locally. A four vector-field with a non- vanishing time component is one 
of the most simple toys that can be used to do so. It will select a preferred direction, a local 
dynamical preferred rest frame at each point of space-time, and can therefore be seen like a 
four-velocity of a fluid present evereywhere (that's why it is called the Aether). Such vectors 
coupling to matter were ruled out by experiments, but one can choose to consider a vector 
that couples only to the metric. A lot of work about such vectors has been done in recent 
years, especially by Kostelecky (|I2]), Jacobson, Mattingly and EUing ([9j, [HI], and [6j ), 
Lim and Carroll ([B] |1]), and Ferreira and Zlosnik f^] [23] I 

2.2 Aether Lagrangian 

Understanding why certain people use a given form of the Lagrangian, and why other use 
something else, include other terms or not requires to look closely at this Lagrangian. What 
follows in this section is an explanation of the logics of the Lagrangians used, after a bib- 
hographic research amongst various papers, especially by WiU and Norvedt (see |2T] [22]). 
Jacobson, Mattingly and Elhng (see [9], [10], and [6]) . 
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To write down an action for a vector-metric theory, one can look for a Lagrangian scalar 
density (scalar so that it is independant of the volume of integration in the action) which 
will, for considerations of simphcity, give a hnear equation for the vector field and second 
order equations at most. The most general Lagrangian scalar density is thus: 

L {A, g) = ao + aiR + a^A'^A^ + agA^A^i? + a^A^'A^R^n + a5(V„A^)(V"A'^) 

+a6 (V^A'^)' + a7(V«A^)(V'^A") + a8(A^V/3A")(A^V-,A,) (11) 

the Oj being constants. All the possible combinations of indices for the quadratic terms in 
covariant derivatives are present here. One can notice that the Aether part of the Lagrangian 
density involves covariant derivatives and thus Christoffel symbols containing first derivatives 
of the metric, so this part of the action will contribute to the metric kinetic terms too. 

But one can simplify this Lagrangian density. First, the difference of the term with the 
coefficient oe and the one with a-j with the one in is Vq (A"V^A^ — V^A") , ie a total 
divergence, which will thus contribute just by a boundary term, according to the Stokes 
theorem. We can therefore choose not to consider the term in 04. 

The exact form of the Aether field depends on what one is looking for. If one wants not 
to keep Lorentz invariance at all or do not mind to do so, one can choose to fix the norm 
of the vector so that it will always have a non- vanishing timelike component and will there- 
fore always be Lorentz-violatin^ The additionnal constraint can be enforced in a maybe 
non-appeahng way, ie using a non-dynamic Lagrange multipher A, but will also simplify the 
equations a lot. 

If one fixes the norm, the terms in 02 and 03 play thus the same role as the one respec- 
tively in ao and a\ and are thus useless. 

The ao term plays simply the role of a cosmological constant (Namely we get a term 
proportionnal to Gap + ^gap by differenciation of the action S = J d'^Xy/^{ao + aiR) if 

A = -|^-)- 

We are left thus left with: 

L {A, g) = a,R + a5(V„A^)(V"A^) + a, [V^A^f + aj{VaAp){V^A") 

+asiAf'\/^A'')iA^V-yAa) (12) 

One can so consider the action: 

S = [ d'xV^ [R + K'^" Vo.A^VpA'' + A (A" A, + 1)] (13) 

K'^" = cig'^^g.a + C25^5f + c,5^y^ + c.A'^A^g,, (14) 

^The Aether will thus carry a non-linear representation of the local Lorentz group because it will not 
belong to a vector-space because its norm being fixed, it will take its value on an hyperboloi'd of the tangent 
space of an event of space-time. 
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This action is the one which was considered by Jacobson, Ehng and Mattingly. Notice 
that dropping the terms in C2 and C4 and considering C3 = — ci, we find K"'^ V aA'^V jsA" = 
^FaaF°"^, where Faa is the antisymmetric Maxwell tensor defined by Faa = VaA^r — V^Aa- 
This simplification was used by Jacobson and by Bekenstein in TeVeS. 

I focused on an action involving a general function F of the Aether kinetic terms. 

S = [ d^x^g [R + M^^ (/C) + A + 1)] (15) 

/C = ^ir"^.V,A^V/3A'^ (16) 

K''^ = c^g'^'^gya + C25°5f + csS^S^^ + c.A^A^g,, (17) 

This action was considered recently by Zlosnik et al. (see [Ml), but with no term in C4 to 
simplify. 

2.3 Fields equations 

I rederived the equations from Zlosnik et al. Aether paper j25j. For I had not studied the 
Lagrangian version of GR before, I had to get used to this formahsm, especially to various 
variations of the action such as w.r.t. the metric and the subtulties of calculations. 

What must be borne in mind when carrying out the variations is that the two dynamical 
degrees of freedom considered are the inverse metric g^" and the contravariant Aether vector 
field A^. The contravariant Aether is chosen (and not the covariant one) just because once 
one has chosen to variate the action w.r.t. g^^ , the result of this variation will be simpler 
seeing the form of K""^ ^„ because we have : 

SA^" , , 

= 18 

whereas: 

9p.p9''' = ^l^ 5^7/.. = -g^.pg.M''' (19) 

and so 1^ = A'^l^ = -g.^A. (20) 
The vector equation is obtained by varying the action w.r.t. A^: 

^ V,(F' J" p) - F'yp = 2\Ap (21) 



5A' 
where: 

. P' =dF 

^ dK 

• is a current: J" a = (i^°^ ^7 +^^" 7a)V/3A^ (One can notice that if iT"^^^ = 
K^"" -ya- like here, for the case considered, J" a = 2K°'^ V/jA''', but defining the current 
this way preserves the generality of the equations for a tensor which does not have such 
a symmetry.) 
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We can get the Lagrange multiplier A from here. 

Variating the action w.r.t. A will give the constraint on the norm: 

A'^Aa = -1 (22) 

For the variation of the action S = J d'^x^—gL w.r.t. the contravariant metric, one must 
notice that = J d'^Xy/^ (sf^ " 25'"/^-^) '^here one uses the fact that: 6g = gg^^^g^^ = 
—ggfifSg^" , g being the determinant of the contravariant metric. 

The symmetry of K""^ simplifies the equations: 



Sg 

with Yai3 a functionnal derivative defined by: 



F'[y^^ + J\-^J^] (23) 



Y^, = V.A^V.A^^^^^ (24) 



The variation of the covariant derivative of the contravariant Aether requires to variate 
the Christoffel symbol (only): 

Sg^P SgC^P SgC^P ^ 

And we have 6 (F^^) = ^ (VaSgpr + 'Vpdgar — '^r^gap) (see Weinberg, [20]) so one even- 
tually find: 

^^^^ = -\v.{:F\J^^'^Ap) - J\^Ap) - J(^p)A")) (26) 

dropping divergence terms which would once more contribute only by boundary terms. 
The brackets denote symmetrization, ie for instance J(a(}) = \ {Ja/s + J/Sa)- 



We have also: 



So putting all of that together, one eventually find: 



- -A^A, (27) 



- a/3 ~ -^0/3 

where T^S*^^^ is the Aether stress-energy tensor: 



G^p = SvrGT"'^"'^" + T^f"'' (28) 



1 



r + ^gapM^:F + XA^Ap (29) 
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2.4 Exploring different regimes 

Once these equations of motions are obtained, one can calculate them for different metrics 
to explore different regimes. 

My supervisor made me consider a Friedan-Robertson-Walker (FRW hereafter) pertubed 
metric such that: 

ds^ = -(1 + 2e(f))dt'^ + a{tf{l + 2eip){dx^ + dy^ + dz^) (30) 

This metric is a perturbed form of the one of a homogeneous and spatially isotropic universe 
(also spatially flat here, ie with no curvature parameter). and are scalar gravitational 
potentials (that are identified in the non-relativistic limit), and a{t) is the cosmic scale factor 
(it can be described as setting the scale of the geometry of space). 

We can recover: 

• the non relativistic hmit by neglecting time derivatives and taking a{t) — 1 to have 

^ _(i ^ 2e0)di2 + (1 + 2ei)){dx^ + dy'^ + dz'^) 

• a homogeneous and isotropic universe by taking e = to have ds^ = —dt^ + a{tY{dx^ + 
dy"^ + dz"^) 

In the following, the equations are developped up to orders in e, but e is not kept for 
a better hghtness. Some equations that are not enlightning to follow the report are put in 
Annex A. 

Einstein tensor Up to hnear order in e we find: 

• Goo = 3i/2 + QHdttl^ - ^Vtl^ 

• Goi = 2{Hdi(P - dtdiij) 

• = (a ' +2a 'a)[-l + 2{(P - iP)] + {dy ^ +d, ^) (0 + V) - Sa^a^ ^ + 2a a dt{(p - 3iP) 

• Gij = -didj{(f) + -0) for i ^ j 

Aether field. Wc take a homogenous and spatially isotropic universe for the background, so 
the Aether must, in the background, respect this isotropy for the modified Einstein equations 
to have solutions, so only the time component can be non zero. The constraint on the norm 
is gapA°'A^ = —1 so in the background, we take A"' — 5q and one can then expand it and 
write: 

A" = (5^ + (31) 
The constraint on the total vector fixes — —(f) with the perturbed form of the metric. 



We can also derive: 
with: 



(32) 
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• Eoi = di<P + dtBi - HE, 

• = [a a +a'^dti> + a d{2ijj — 0)] Sij + diBj 
The other components vanish. 

Matter. For matter fields, we can take: 

j^rnatter ^ ^ p)^^^^ ^ p^^^ (33) 

which is the stress-tensor of a perfect fluid without any anisotropic stress, with a density p, 
a pressure P and with the fluid four- velocity satisfying g^uU^u^ = — 1- If we consider a non- 
relativistic fluid, ie with no spatial components of u^, for this metric: u^j, = (—1 — e0, 0, 0, 0). 

With this metric we have thus: 

To7«-^- = (l + 20)p (34) 

rpmatter rpmatter q (35) 

7;7*te^ = a2(l + 27/;)P (36) 

I carried out the calculations of the Einstein equations up to linear order analytically 
with this metric, except for the cross terms, and for the calculations are very tedious, time- 
consuming, and it is really easy to make mistakes, I tried to compute them with a calculus 
software, to check them and to calculate the cross-terms. I tried lots of Mathematica packages 
("tensorial", "GREAT", "GRTesting", "xtensor") before flnding the Maple "tensor" package 
that is a free one which can really carry out calculations with tensors, especially in which one 
can really enter components of a tensor (the Aether here) and that can compute covariant 
derivatives of it. I thus wrote the Maple sheet that is in the Annex B at the end of this 
report. 

2.4.1 Static limit 

In the static hmit, the spatial terms of the Aether appear only at second order in all the 
equations. 

T^*'*'^'' has no cross-terms (up to hnear order), so we find: 

G,. = -9,9,(0 + ^1;) fort^j I ^ = -V^ (37) 

rp Aether _ ipmatter _ q Jqj' i ^ j J 

So the only non-zero component of the Einstein tensor is: 

Goo = 2V'0 (38) 
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We have: 

T^f'''' = (C4 - ci) V. (F'V0) - ^ (1 + 2£0) FM^ (39) 
so we find the equation 

V. ((2 + (ci-C4)F')V0) = 87rG'p (40) 
It looks tempting to recover MOND with such a Poisson equation including an extra term. 

The kinetic scalar K is in this case, up to second order: 



Restricting to the case of C4 = considered by Zlosnik et al. ([H]), we see that we can 
recover a MONDian regime in which V. (|V0| V0) oc p in the limit of small |V0| writing 
lim [2 + (ci) F'] oc K2 so we can get a real MONDian limit if we consider M of the same 

|V0|«A/^ 

ordre of Oq. |^ More details and discussions are given in [21]. 

2.4,2 Homogenous and isotropic universe 

The only non-zero component of the Einstein tensor are: 

• Goo = 3i/^ 

• Gxx = Gyy = Gzz = — (o^+2aci) 
We find: 

jAether ^ 3^^ (2F' + F' H + F'^ (42) 

The vector equations gives the Lagrange multiplier: 

A = 3 (ci + C2 + C3) F'- - 3c2 [F' H + F'-) (43) 

SO that we get: 

rj^Aether ^ 3aF'^ - ^M'^ F (44) 
where a is defined as a = ci + 3c2 + C3. 
We have therefore the 00 modified Einstein equation: 

3{l -aF')H^ + ^M^F = 87iGp (45) 



^This was done by Zlosnik et al. (see |2l]). With no C4, K is ensured to be positive for consistence of 
a quantisized theory of the Aether carried out by Lim in [13]. Including the C4 safely would require further 
examination. 
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We can also calculate: 

r^Aether ^ j^Aether ^ j^Aether ^ (^a+aO^-aF'aa + ^FM^ (46) 

and so we get the modified pressure equation: 

- {I -2aF')H^ -2(1- -aF'\ - + a F' H - -FM^ = 8nGP (47) 

\ ^ J Qj 2 

Zlosnik et al. identified the additional terms that appear here with a cosmological con- 
stant. More details and discussions are once again present in 
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2.4.3 General perturbed metric 

Aether with no spatial terms We consider first the case where the only non-vanishing 
component of A is Aq, the spatial components being zero. 

Kinetic scalar We get for K: 

SaH'^ 6aH jj, , ^ ,^ , 1 (c4 - Ci) 2 

+ ^ {-AHi^dti^ + 5i/V' - ^H<t>dti^ + {dti^f) + (48) 

Density and pressure equations The Too component of the Aether stress-energy 
tensor is: 



rj^Mther ^ (C3 + C4 C,) ^ 



+3C2 



F'{2H'^ + - + QHdtil^ - Hdt(t> + 9t 2 V) + F'iH + 9^^) 

/7 



- ^(1 + 2(/))M2f + AAoAo 



(49) 

The vector equation gives: 
A = 3 (ci + C2 + C3) F'//2 _ I^F'^ + F'H^- {F'di(t>) + 3 (2 (ci + C3) F'i/ - F') 9*^' 

+3c2F'Hdt(j) - 3c2F'd^ij + 6 (^C2 (^F'^ + F' i/^ + (ci + C2 + C3) F'H^^ 

(50) 

Hence we have: 

j^Aether ^ _^l^v.(F'V0) + 3aF' + QaF'Hdt^ - ^(1 + 20)714"'^ (51) 
and thus: 

3(l-Q;F')ii'^-^VV+^^^^V.(F'V0) + 6(l - aF') + -(1 + 2</))M2F = STrGTo™""^" 

(52) 

We have T^"-^^'^'^ = (1 + 20) p, so we get the modified density equation: 

3(l-Q;F')i/^-4vV+^^^^^V.(F'V(?))+6 (I - aF') Hdtip-6 (1 - aF') H^(f)+^M^F = SttGp 

(53) 
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The spatial diagonal terms are equal and we have for instance: 

+a -3 (2F'a a + F'a^) + F'a a dt<p - F'a^d^ip + 2 (^F' (2a +aaj + F' a aj (0 - ^) 

(54) 

We find therefore, as for matter, TT^""^'' = (1 + 2ip) P, the modified pressure equation: 

- (1 - 2aF') H'^-2[l- laF') - + aF'H- \m'^F + (0 + -0) 

+ (-6H (1 - aF') + a F'^ dti/j + (2 - aF') Hdtcl) - (2 - aF') 9^^^ 

+ (^-2a F'H + 4{1- aF') ^ + 2(1- 2aF') H^^ = SnGP (55) 

Cross terms We have 

Tot ^'^^ = (C4 - ci) (F'a^a.^ + (f' +/7F') (56) 

Thus we have, since the T^'***^'' is taken to be zero here (for a non-relativistic perfect 
fluid): 

(-2 + (ci - C4) F') ^t^i^|J + [(2 + (ci - C4) F') + (ci - C4) F'] 9^0 = (57) 
We have for i ^ j, up to hnear order 

j^Aether ^ g (58) 

So since the stress-tensor of matter has no cross terms either in the case considered, given 
the form of the Einstein-tensor cross terms Gij = —didj{4> + ^) for i j, we can, as in the 
static case, identify the potentials cf) and -0. 



Aether with spatial terms We consider now also spatial terms: we take more precisely 
a covariant Aether of the form (—1 — e0, eBj., eBy, eBz). 

Kinetic scalar We have, up to linear order: 

The expression of K to second order is very long and not very useful such. It is in Annex 

A. 
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Density and pressure equations We can derive: 

rpAether ^ - ~f\ v .{F'V(j))+diF' dtBi)+3aF' H^+6aF' HdtiJ+2a^F' diBi--{l+2(l))M^ F 
a'' o'^ 2 

(60) 

The expression of the Lagrange multiplier A used to get this is in Annex A. 
We have thus: 

3(1 - aF')H^ - ^V^i/j + ^^^(V.(F'V0) + d.F'dtBA + 6(1 - aF')Hdtij 



-2a-^F'diBi + J(l + 2(t>)M''F = SttGToo (61) 



and: 



3(1 - aF')H^ - + ^^^(V.(F'V0) + d^F' .dtB,) - 6(1 - aF')H\ 



-,F'd,B, + \ 
2 

We find: 



+6(1 - aF')Hdtij - 2a—F'diBi + -M^F = SirGp (62) 



T-- = .(-r(2^.aa)-r„a, + la^(1.2.,™^ 

-3 (2F'a a + F'a^^ dti/j + F'a a dtcf) - F'a'^d^x/j + 2 (^F' (2a +ad^ +F'a aj (0 - x/j) 

- (ci + C3) (f' +F' a) d^B^ - (c2 F' + (ci + 4c2 + C3) F' a) 9,5, - a a 

- (ci + C3) F'^i^.E, - caF'^t^i^, 

(63) 

The pressure equation becomes: 

- (1 - 2q;F') - 2 (^1 - ^aF'^ '^ + aF'H- ^M^F + ^ v' (0 + ^) 
+ (^-6// (1 - aF') + a F') + (2 - aF') Hdt(t) - (2 - aF') d'^i^ 



+ ( -2aF'i/ + +4(l-aF')- + 2(l-2aF')i/^ 



a / F' +AF' a 1 1 

+ ^ 3 ^^^^ + " ^^^'^^ + g^'^t^i^i I = 87tGP (64) 



Cross terms We find: 

To, = (C4 - ci) (F'9t9,(/. + (f' +i7F') 9,0) 
(C4 - ci) (F'9t 2 5, + (F' +F'H)dtB,dtB,^ + a (^F'^ + F' - F'i/^^ 



+ ^^^ (5^F' {d^Bi - d,B,) + F'a, {d,B, - OA)) (65) 
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and so we have: 



(C4-C1) (^F'dt^B, 



(-2 + (ci 



X 




(66) 



We find that the Aether stress-energy tensor has now spatial cross terms in linear order. 
For z 7^ j, we have: 



We see thus that the spatial terms of the Aether makes the two scalar potentials differ 
from one another. We cannot therefore identify them as we could do for an Aether without 
spatial terms. Notice that in GR, these two potentials can only be made to differ by off- 
diagonal terms of the matter stress-energy tensor, ie by anisotropic stress. One can also see 
that Tij is zero if Ci = — C3, in the magnetic case. 

I have thus got the perturbations of the Aether stress-energy tensor and the Einstein 
equation for an Aether with or without spatial terms and a Lagrangian involving a general 
function of the kinetic term K. 

These equations are very general. Simplifying assumptions must be made for physics to 
arise. I unfortunately did not have time to go on much with this during my internship but 
work is currently being done about it, keeping this general perturbed metric, but considering 
for instance the vector perturbations as gradient of a scalar field. For instance, in this 
particular case, and considering F' as constant, the the Ox equation, would reduce to the one 
of a harmonic oscillator, with a negative string contant and a source term which is a function 
of the scalar potentials and ip. One can therefore explore (numerically) the growing of this 
scalar potential from an initial instability. Similar work was done by Dodelson and Liguori 
(see [5]) , but for the Bekenstein TeVeS Lagrangian. The generality of the equations presented 
here gives nevertheless the opportunity of exploring various cases just by simplifying them. 



lAether 



Cl + C3 

2 



{F'Hd^.Bj) + F'dtd^^Bj) + F' 9(,5,)) 



(67) 
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2.5 Scalar field 

We can try to take F{K) = K and make the coefficients of the Lagrangian depend on a 
scalar field "instead". Similar work has been done by Kanno and Soda ([H], [11]) in the 
context of cosmic infiation. 



With: 

S = I d'^x^/^ ' 



WtiG 



+ L{A,g,x) 



N 



+ s 



M 



and: 

L{A, g, x) = tttV f + -^(^"^^ + l)-Vix)- ] (69) 

the scalar equation is: 

_dy^ ^ V.A^V"A, + ^ {V^A-f + '^V^A^V.A- + '^A-A^V^A.VpA^ = -V.V^x 
dx dx dx dx dx 

(70) 

and the Einstein equation, considering (7"^, A™" and the scalar field as three degrees of 
freedom is: 

I defined here (in the action) the scalar tensor by analogy with the Aether tensor as far 
as the coefficients in the Einstein equations are concerned. 
We have: 

^.scalar ^ (^(^ ^ y ^^^^ ^ ^V.^V^X (72) 

Cosmology We consider the case of a homogeneous and isotropic universe with: ds'^ = 
—dt^ + a{tY{dx^ + dy"^ + dz^). The scalar field can depend only on time. (If we make 
it depend on spatial coordinates, the scalar stress-energy tensor will have non-cross terms 
rpscalar ^ 1 (^.-^^ (djx) for i j that are zero using the Einstein equations.) 

The only non-vanishing terms of the scalar stress-energy tensor are thus the diagonal 
ones, which are: 

T^t''' = \id,x? + lvix) (73) 

T:f^' = T^f^^ = T:f'^^ = --y (- {o^xf + 2v {x)) (74) 



The scalar equation is: 



dV nda „o „ / s 

+ 3H'—=dh + Smx 75 
dx dx 



with a (x) = ci (x) + 3c2 (x) + Ci (x)- 
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and the Einstein equations are: 

3 (l - ^) = i {d,xf + \V ix) + SttGp (76) 

- (l - ^) - 2 (l - (X)) l = \ {daf - \V (X) + SvrGP (77) 

If we try to take a constant scalar field and a simple scalar potential V (x) = \M'^x^-, the 
scalar equation is so: 

-M\ + 3i/2^=0 (78) 
dx 



and the Einstein equations are: 



3(1- ^ 1 = -M'x' + 8nGp (79) 



'i-^)i/^-2^-l«M^^^ = -i 

or, using the scalar equation: 



2 2"(x)l^ = -^MV + 87rGP (80) 



-^).-.(:-i"w)MS(|)^-'- 

We can try to identify the additionnal terms in the modified Einstein equations with a 
cosmological constant A (hke what can be done for a simple Aether, with no scalar field [24j) 
ie such that Gajs + ^Qap = SttCTq/j, like what was done in Zlosnik. We see that we should 
have: 

We can't therefore do it for H that is time-dependent appears in the coefficients of the 
differential equation whereas a is a function of the sole scalar field which is taken constant 
here. 
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2.6 Maple code 

The Maple sheet I wrote aims at computing the equations of the report with a large choice 
of parameters. It follows the exact notations (indices, etc) and the exact steps of the report. 
Maple can't carry out the calculations from the Lagrangian, it can't variate the action w.r.t. 
the metric or the Aether, so the tensorial equations after differenciation must be entered in it. 

The sheet of the Annex is expanded for a simple FRW metric because it is the metric 
that gives the lightest Maple results. These results can, in spite of collection of terms by 
Maple, be several pages long for a FRW perturbed metric and an Aether with spatial terms. 

The first section allows to choose a system of coordinates and a metric. Five types of 
metric are displayed in the annex: a simple FRW metric, a perturbed metric with static 
potentials, the FRW perturbed metric of the report, the same in conformal time (conformal 
time r is sometimes used in papers, it is defined by airYdr'^ = dt^, and a spherically symetric 
metric which I added to the Maple code but did not mention in my report. One can add 
off-diagonal terms easily in a metric, for all components are defined separatly. All tools to 
compute the Einstein tensor, and the Einstein tensor are computed (Those are procedures 
included in the package). The Aether field is then defined, including spatial terms or not 
(when it is consistent with the chosen metric, ie not for a FRW metric for which spatial terms 
will break the isotropy). 

The second session aims at computing the Aether stress-energy tensor. All intermediate 
tensors are defined so that one can change them individually. The tensor K^^ can be easily 
changed if one wants for instance to put no constraint on the norm of the Aether and thus 
to add terms in it. One can choose to take a general function F of the kinetic term K, to 
take F{K) = K, to have a "simple" Aether theory, to take a Maxwell-like K to simplify even 
more, and one could enter a specific function. Some components are expanded at the end of 
the section. Notice that the input lines are very long because the results given by Maple are 
often messy, so one must arrange them by collecting terms (in the Annex, some commands 
used for other metrics are left in the input lines). 

The third section computes the stress-energy tensor (the one of a perfect fiuid) for matter, 
given the four-velocity defined in the first section. 

The fourth section is optionnal. I wrote it to include a scalar field whom the coefficients 
of the Aether Lagrangian depend on (see last section). It computes the scalar stress-energy 
tensor and the scalar equation. 

The fifth section computes the modified Einstein equations. 

This calculus sheet was approved by Tom Zlosnik that works on the Aether in the Uni- 
versity of Oxford. (I sent it to him but did not collaborate further with him.) 
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3 The EFE: example of a galaxy disc with a bulge 

EFE: External Field Effect 

The MOND gravitationnal equation is such that if any external field is considered, it will 
break the Strong Equivalence Principle which reads that gravity can always be replaced by 
an accelerating frame. This is fundamental for GR, but once again, it is just a principle and 
might not be true. Nature might not be as simple as we often wish it is. The external field 
in MOND has effects on the interior dynamics of a system embedded in it. This should not 
be confused with tidal effects because this occurs even if the external field is uniform. 

I studied an example of how a static external field can affect a system, a galaxy disk with 
a bulge around it. This real MOND problem involves various galactic dynamics subjects: 
potential issues, phase space densities, etc... My supervisor told me about his work on 
this problem: I TeXed it for him, as well as another two pages about phase space density 
considerations, explaining it with clarified notations, drawings, etc... 

3.1 Kuzmin disks and the Plummer model in Newtonian gravity 

The Plummer model was first used by Plummer to describe globular clusters, massive stellar 
systems containing 10"^ — 10^ stars in a nearly spherical distribution (see [3]). 
We consider the spherical potential: 

GM 

^Plummer — , (84) 

rc is the core radius (the radius where the surface brightness has fallen to half its central 
value). 



We find thus: 



d f 2d^P 



PP/_(r) = ^,^{r—] (85) 



3Mr 



2 



Att (r2 + r2) 2 



(86) 



3Mr 



2 



47r(GM)'*''^"~ ^^^^ 

To describe galaxy infinitely thin disks, Kuzmin introduced a potential pair, consisting 
of two fictive point masses. We consider a disk and an axisymmetric potential given by: 

GM 



K 



We have therefore: 



GM 



GM 



R2 + (\z\ + z,y 

above the plane 
below the pla 
on the plane 



= ) — ^ below the plane (^gg) 
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^ - (Ft.z) 



2^ .. 



galaxy disc 



Figure 1: Kuzmin potential 

So above the plane, the potential is the one of a (fictive) point mass located at —zq, and 
below the plane, the potential is the one created by a fictive point mass at zq- 

The force has therefore a z component discontinuous on the plane of the disk (the R one 
is continuous). The density is zero above and below the plane, because the potential is there 
"point-mass hke", but one can apply the Gauss theorem to an infinitely thin box for instance 
with upper and lower parts right above and below the plane to calculate the surface density 
a of the disk: 



9$ 



K 



dz 



U=o+" 
1 



K 



dz 
Mzo 



2=0- 



(i?2 + ^2)2 



(90) 
(91) 



We have thus a mass-distribution which is purely in the disk. 

If we want now to consider not only a thin disk, but a bulge as well, we can take now the 
potential: 

GM , , 

, (92) 

^B? + Tl + {\z\ + z^f 

so that we get spherical Plummer potentials above and below the plane: 



GM 



GM 



above the plane 
below the plane 



(93) 



so we have a Plummer density above and below the plane. 



We can once again use the discontinuity of the potential at 2; = to calculate the surface 
density of the disk, which depends now on R: 



E{R) 



Mzo 



(94) 
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is: 



Its mass is: 

POO 

Md^sk = / 2TTj:{R)RdR (95) 
Jo 

Mzq , , 

and with p+ the density above the disk, because of the symmetry, the mass of the bulge 

/•OO POO 

Mkuige = 2 / 27rp+ {R, z) RdRdz (97) 

3.2 Kuzmin disk with a bulge in MOND 
3.2.1 With no external field 

To carry out calculations, we assume the Bekenstein /i function ( [I]) : 

/i (x) = ^ with X = — (99) 

The MOND equation for a mere sphere is: 

V.{^^{x)V^)=A'KGp = V^^N (100) 
^ V. (/i(x)V$- V$^) = (101) 

In general, this thus implies: 

/i(x)g = + V X h (102) 

But for spherical symmetry, we have (one can see it by integrating and using the Stokes 
theorem): 

/i(x)g = gTv (103) 
With the particular form of the /i function, we get: 

\g\ = \gN\ + a/Moo (104) 

We can consider a gravity center at the location (0, — -Zo), below the plane, and we have 
thus: 

- \gN (r+) I e,.+ - a/ Qn {r+) aoe^+ (105) 



+) 



The gravitationnal field Q]^ (r+) is a Newtonian one, so a Kuzmin one here, ie: 

d ( GM , , , 

g^M = --^\ = (106) 



dr 

GM [R + (z + zo 



R^ + r^, + i\z\ + zof 



{R^ + r^, + {\z\+zoy) 



2x3/2 



(107) 
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and the gravitationnal field below the plane is the mirror image of this. 
The potential itself is: 

/•r+ 

= - / 9{r+)dr+ 
J +00 

[gN (r+) - y/\gN{r+)\ ao) dr+ 



' +00 

— ^Newton + ^effective DM 

3.2.2 Sphere dominated by an external field 

We consider now a sphere dominated by an external field in a MOND gravity, ie : 



(108) 

(109) 
(110) 



V. ( ^ ( ^ 1 V$ 



with - V$ = gext - V$i„t , Igextl > |V$ 



int 



(111) 

(112) 



We can see that the "external field effect" of MOND is due to the non-hnearity of this 
equation: we can expand /i: 



I gext - V$ 



int 



ao 

ext I \ dp 

dg 



ao 



gextV$int 



Igextl 



int 



Igextl dZ 

and the MOND equation is thus: 



with He = P 



I gext I 

ao 



and Lf. 



din jjL 



din |g| 



(113) 
(114) 

(115) 



Igextl 



V. 



Igextl dZ 



;-gext + V$int) 



(116) 



so up to hnear order in ' , ^ , ' , we have. 

gext 



-V.gext + V^^int + L, 



int 



The first term is zero for the external field is taken uniform, so: 



(117) 



92 



Pe 



+ 



92 



|2 H 



^dz^ 



int 



AnGp 



dZ^ 



$int = A^t:G—, with Li = Le + l 



(118) 
(119) 
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For a Newtonian gravity (with fi = 1 and Lg = 0) we would find the mere Poisson equation 
for the internal potential, whereas here, we find an equation similar to the Poisson equation, 
but anisotropic and for a higher density ^. 

With the transformation 

X' = X, y = y, z' = —= (120) 

V Li 

we find the Poisson equation in these new coordinates, still with the same higher density: 

Vf^int = AttG-^ (121) 

where the , refers to the {X',Y',Z') coordinates. 

3.2.3 Ellipsoid dominated by an external field 

We can consider the physical density: 

P = -T=PPlummer{r') (122) 
V -^1 
1 "^Mr^ 

" , with r' = + r'2 + (123) 



'Li (r'2 + r2) 2 



where the factor is such as to get a total mass integrated over all space equal to 
M. This is the density of an ellipsoid (the surfaces of equidensity are surfaces such that 



X^ + Y'^ + = cste, ie elhpsoids). We have the perturbed MOND equation 



V'^int = 47rG^ (124) 



/i 



e 



= 471^ ^"" ' (125) 



(r') = ^'''^ (126) 



So we have: 



The relation between the density p and the potential (pint differs therefore by a factor 

,,5r5/2. 

5 SMr2 

p = -pilL! \ (127) 
' 471(^71/)^ ^ ^ 

If we assume that the origins of the elliptic potentials are in the {X' , Y') plane, at 
(— Xq, — Fq') and (Xq,Fq'), that the external field is also in this plane, and we define the 
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coordinates {x, y, z) such that the {x, z) is the same as the {X, Z) one and by requiring that 
the potential is continuous on a plane 2; = 0, we have on this plane: 



CM 



1 



int 



r' = r' 



GM 1 



so that wc can write for a point on this plane: 

{X' -x'^f + V' + ^z' -z'^f = 



(128) 
(129) 



(130) 
(131) 




disc 



Figure 2: Angles and coordinates 

With the angle between the z axis and the Z one (ie the one along wich the external 
field is), we have therefore: 



tan df. 



Xn 



If we define /3o and ^0 such that: 



tan /So 
and $0 



Li- 



Xo 



(132) 
(133) 



(134) 
(135) 
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we find: 



tan^e - tan/3o , . 

tant^o = tt (136) 

_ (Li - l)sin6'e cos6'e . . 

,2/1 U'^'J 



Li cos^ 6*6 + sin^ 9, 



e 



So as we turn on an external field, we see that the center of the potentials shift and their 
axis rotates. 

We have above the plane: 

1 GM , , 

^int = , = (138) 

y Lir2 + A (x + x^f + Liy2 + S (^ + z^f + C (x + Xq) (^ + ^o) 

with ^ = (Li cos^ + sin^ ^e) , (139) 

S = (Li sin^ + cos^ ^e) , (140) 

C = 2(Li - l)coseeSin^e (141) 



and we can also find 

l + Lcos2 



L cos 6'e sin 6le 

^0 = -^OT— -r— ^ (142) 



3.2.4 Potential energy 

I calculated the potential energy for the potential: 



PE^^ 1 1 1 p<j),^,d^x (143) 




We find: 



3GM^ / TT zoVA 1 ^0 , . 

PE = — ^ arctan — 144 



>1 = Li cos^ + sin^ (145) 

This sole result is not enlightning. With no external field, we would have just the first 
term. One could explore Virial theorem related questions, carry out some numerical simula- 
tions, etc. I did not go much further about this because I came back to the Aether. 
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Conclusion 

During this internship, I discovered MOND, its birth, building, and got in touch with some 
MOND issues. I also discovered the Aether, which is a field by itself and is finked to a 
fiuge amount of various attempts of modifications of General Relativity. I rederived the 
fundamental equations of the theory, calculated them for different metrics and found a right 
way to compute them. I calculated them for a FRW perturbed metric, and a general function 
of the kinetic terms appearing in the Lagrangian, which had not been done such yet. I wrote 
a Maple code (approved by Tom Zlosnik of the University of Oxford) that can be used for 
various computations around the modified Einstein equations. This sheet is a useful tool 
to carry out calculations and can be easily modified. I also got interested in a problem of 
galactic dynamics in MOND, involving the special "external field effect". Some further work 
is currently being done about the Aether and the growth of perturbations. 




^The pictures of the title page and the conclusion are pictures of WestSands, a beach of St- Andrews (the 
town were I made my internship). They have no special meaning, they are just put here for their beauty. 
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